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ABSTRACT 

We show that near-infrared observations of the red side of the Lya line from a 
single gamma ray burst (GRB) afterglow cannot be used to constrain the global neutral 
fraction of the intergalactic medium (IGM), xjj, at the GRB's redshift to better than 
5xh ~ 0.3. Some GRB sight-lines will encounter more neutral hydrogen than others at 
fixed xh owing to the patchiness of reionisation. GRBs during the epoch of reionisation 
will often bear no discernible signature of a neutral IGM in their afterglow spectra. We 
discuss the constraints on xh from the z = 6.3 burst, GRB050904, and quantify the 
probability of detecting a neutral IGM using future spectroscopic observations of high- 
rcdshift, near-infrared GRB afterglows. Assuming an observation with signal-to- noise 
similar to the Subaru FOCAS spectrum of GRB050904 and that the column density 
distribution of damped Lya absorbers is the same as measured at lower redshifts, a 
GRB from an epoch when xjj = 0.5 can be used to detect a partly neutral IGM at 
97% confidence level « 10% of the time (and, for an observation with three times the 
sensitivity, 30% of the time). 
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1 INTRODUCTION 

Tomorrow, a gamma ray burst (GRB) may be observed that 
originates from the death of one of the first stars, during 
the epoch of reionisation. Despite the great distance to this 
burst, it will be the brightest gamma ray source on the sky 
for several tens of seconds, one of the brightest cosmological 
X-ray sources for hours, and its afterglow will be observable 
for weeks in the near-infrared (and, for the first few hours, 
brighter than any z ~ 6 QSO). Much of the optical and 
near-infrared light will be obscured by the Lycv forest, and 
this obscuration will enable the strongest constraint to date 
on the neutral hydrogen fraction of the intergalactic medium 
(IGM) at the burst's redshift. 

In fact, such an occurrence may alread y have 



been realised. lHaislip et all il2006f). iKawai et all (l2006h . 
iTagliaferri et all (|2005f )~ and iTotani et all (|200rj ) ob- 



served and analysed the optical/near-infrared afterglow of 
GRB050904, identified to be at z — 6.3 - possibly during 
the reionisation epoch and the GR B with the highest con- 
firmed redshift. ITotani et all (|2006l ) derived the constraint 
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on the global neutral fraction xh < 0.6 at z = 6.3. In this 
paper, we discuss the assumptions that went into their anal- 
ysis, and we investigate how realistic modelling of reionisa- 
tion can affect constraints on xh from GRB050904 and from 
future 2 > 6 GRBs. 

The Swift satellite has greatly increased the sam- 
ple of GRBs with known redshifts in the last two 
years dGehrels et all I2004D. Future missions such as 
EXIST dGrindlav fc The Exist Team] l200rj) and JWST 
^Gardner et all l2006h will further enhance our ability 
to detect high-redshift GRBs and will enable more de- 
tailed follow-up studies of their near-infrared afterglows. 
Interestingly, approximately one-half of Swift bursts are 
"dark bursts" - bursts that have detected X-ray after- 
glows, but that hav e no measurable optical emission (e.g., 
iFilliatre et all |2006l) . While it is probable that most dark 
bursts originate from low-redshift, dust-rich galaxies, a frac- 
tion of dark bursts may originate from z > 6 and are "dark" 
because Lya absorption fr om the high-redshift IGM absorbs 
the optical emission fe.g.. lMalesani et alj|2005h . 

In addition to their extreme luminosity, there are several 
other advantages to studying reionisation with GRBs com- 
pared to other probes of this epoch. First, the afterglows of 
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high-redshift GRBs are observed at earlier (brighter) times 
in the source frame than those at lower redshifts, so the dim- 
ming owing to increased luminosity distance is nearly can- 
celled , and the observed flux is almost independe nt of red- 
shift (ILamb fc ReichartJlioOll ; ICiardi fc LoetfeOOOl ) . Second, 
unlike the spectra of galaxies and quasars, the intrinsic after- 
glow spectrum of a GRB is a featureless power-law at the rel- 
evant wavelengths, allowing a more precise measurement o f 
absorption owing to a neutral IGM ijBarkana fc L ocb 2004). 
Finally, since the theoretical expectation is that most of the 
star formation at z > 6 occurs in halos with m ~ 10 9 Mq 
and because observations at z > 6 currently probe only the 
most massive galaxies and QSOs (m > 10 11 Mq), high- 
redshift GRB host galaxies should be less massive than 
galaxies selected in another manner. Consequently, GRB 
host galaxies will sit in smaller HII regions during reionisa- 
tion (on average) than galaxies selected by different means. 
Therefore, GRBs will suffer a larger Lya IGM absorption 
feature. 

In this work, we do not concentrate on wavelengths 
blueward of source- frame Lya (in the Lya forest) to de- 
rive constraints from GRBs. Any blueward flux indicates 
the presence of ionised gas at the redshift of the transmis- 
sion. However, at high redshifts there is little or no flux in 
the Lya forest, even in ionised regions, owing to the increase 
in density with increasing redshift, the decrease in the size 
and in the number of voids, and the decrease in the am pli- 
tude of the ionis ing background (e.g-. lBecker et al1l2007l and 
iLidz et all [20071 ). As a result, it is difficult to distinguish a 
partly ionised IGM from a fully ionised one with the z > 6 
forest (jBecker et al.ll2007l ; lEidz et all 120071 ). Future obser- 
vations of the z > 6 Lya forest from additional QSOs and 
GRBs will aid reionisation studies, but is unclear whether 
such studies will ever provide definitive evidence for neutral 
pockets in the IGM. In contrast, the shape of the line pro- 
file redward of Lya is sensitive to a substantially neutral 
IGM and, t herefore, can be used to unambiguously detect 
reionisation l|Miralda-Escude|[l998l ) . 

Little is known about the rate of GRBs at z > 6. We 
assume that the rate of long GRBs traces the massive star 
formation rate (SFR) for most calculations in this workQ 
The assumption that the GRB rate traces the massive SFR 
is suppo rted by observations of lower redshift GR B host 
galax ies llBloom et alj |2002| ; iDiorgovski et all l200ll ) . How- 
ever, iKistler et al.l (|2008l ) found that the GRB rate is four 
times higher at z « 4 than if the GRB rate exactly traces 
the SFR. Other properties of a galaxy apart from its mas- 
sive SFR might be co r relate d with its rate of GRBs. For 
example, IStanek et ail (|2006h found that z < 0.25 GRBs - 
GRBs that are typically under-luminous - are preferentially 
in metal-poor galaxies. 

Making the a ssumption that t he G RB rate traces 
the observed SFR, Sal vaterra et alj ((2008) predicted that 
SWIFT will be triggered by 1 — 4 bursts a year above 
z = fQ and that the EXIST mission would observe 10 — 60 

1 We do not consider the other class of GRBs, the "short" GRBs, 
in this study. These bursts, while still cosmological, are more local 
than long GRBs and typically do not have a detected afterglow. 

2 This rate is larger than the SWIFT rate of 1 identified z > 6 
bursts in three years, but it is possible that these numbers can be 
reconciled in light of these dark bursts. 



bursts a year. Other studies have predicted even larger rates 
l|Bromm fc Loebll2002l ; iDaigne et aLll2006l ). 

At z > 6, POPIII stars with average masses of ~ 
100 Mq may exist. It is unclear wh ether the d e ath o f 
a POPIII star can result in a GRB. iFrver et all (|200ll ) 
identified a mechanism that might produce GRBs from 
POPIII stars. However, once the interstellar metallicity 
reaches a critical value of ~ 10 -3 ' 5 solar in a high-redshift 
galaxy, POPIII star formation quenches an d the normal 
mode of POPH star form ation begins (e.g., iMackev et al.l 
l2003t lYoshida et alJliooi ). and this mode is known to pro- 
duce GRBs. Most if not all of reioni sation likely owes to 
photons from PO PII-like stars (e.g., [Sokasian et al. 1 12004 
iTrac fc Cenll2007l ). 

In Section[2] we discuss the absorption profiles of a neu- 
tral IGM as well as of a damped Lya absorber (DLA), and, 
in Section [3] we discuss our simulations of reionisation and 
their implications for the amount of IGM absorption in a 
GRB afterglow spectrum. Section|4]quantifies the detectabil- 
ity of a neutral IGM from GRB afterglow spectra, and Sec- 
tion [5] describes the constraints on xh from the z — 6.3 
burst, GRB050904. For our calculations, we adopt a cos- 
mology with Q m = 0.27, Ov = 0.73, n b = 0.46, a 8 = 0.8, 
n — 1, and h — 0.7, which is consistent with the most recent 
cosmic microwave ba ckground and large scale structure data 
l|Spergel et alj|2007l) . We express all distances in comoving 
units unless otherwise noted. 

When this proj ect was nearing comple t ion, w e learned 
of a similar effort bv lMesinger fc Furlanettol (|2007l ) and refer 
the reader there for a complementary discussion. 



2 THE RED DAMPING WING 

In the standard picture, GRB afterglows result from 
shells of relativistic matter colliding with the interstel- 
lar m edium, shocking , and radiating via synchrotron emis- 
sion (|Sari et all Il998l ). The observed flux of a GRB just 
redward of GRB-frame Lya is given roughly by T — 
A (v/va)^ 13 exp[— Ta.{v)\, where A is a time-dependent am- 
plitude, f3 is a time-dependent power-law index, v a is the 
frequency at the Lya line centre, and 

. . f Zg dz c , . / v(z).\ 

+ iV HI a a (, ZDLA (l + ^l)). (1) 

Here, v z = v (1 + z), z g [zdla] is the redshift of the GRB 
host galaxy [DLA], a a (y) is the Lorentzian-like Lya natu- 
rally broadened scattering cross section at frequency u, and 
v(z) is the line-of-sight peculiar velocity of the gas. The first 
term on the right-hand side of Equation (JlJ owes to IGM ab- 
sorption and the second owes to absorption by a DLA with 
column density Nm- The extended absorption profile red- 
ward of Lya shaped by the factor exp[— T a (y)] is commonly 
referred to as the "damping wing" . 

Let us first ignore absorption owing to a DLA and the 
effects of peculiar velocities. For an isolated bubble of co- 
moving size Rb in a homogeneous medium with neutral 
fraction xh, Equation (QJ simplifies to the expression (e.g., 
lLoeb fc Rvbickll 19991 ) 
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Figure 1. Illustration of the GRB afterglow flux near the GRB- 
frame Lycr wavelength. The top panel depicts the afterglow flux 
for both an ionised universe (dashed curves) and with xh =0.5 
and Rb = 10 Mpc (solid curves). The bottom panel illustrates the 
effect of bubble size on the afterglow flux, assuming xjj = 0.5 and 
iVei = 10 20 cm -2 . The normalisation of the flux in these panels 
is arbitrary, and it has been adjusted in the top panel to separate 
the sets of curves. These curves are calculated with /3 = 1.25 and 
z g = 6.3, and A is the observed wavelength. 




Figure 2. Slices through the middle of the simulation box for the 
three models discussed in the text, each with width 186 Mpc and 
with 6.9 < z < 8.5. The white regions are fully ionised and the 
black are fully neutral. The top row uses snapshots with volume- 
weighted neutral fractions of xh = 0.7, the middle uses those 
with xh = 0.5, and the bottom uses those with xh = 0.2. 



T a {v) « 900 kms" 1 x x H f ^ 9 J 

x (EM^±_ c !±^y\ (2) 

V (1 + %) f« J 

or r(i/ a /(l + z g )) ~ xh for a 1 proper Mpc bubble, noting 
that H(z g = 7) « 900 km s _1 proper Mpc _1 The width 
of the damping wing feature from a neutral IGM (defined 
here as the the wavelength for which exp[— T a (v)] < 0.9) 
is ~ 25 (1 + Zg) A. The width of this feature is typically 
broader than the feature due to a DLA, which falls off as 
A\ a ~ 2 rather than the AAq. -1 scaling for IGM absorption, 
where AX a is the difference bet ween a given wavelengt h and 
that of Lya in the GRB-frame (|Miralda-EscudeHl998r i . 

Figure [T] illustrates the effect of the column density of 
the DLA (top panel) and the size of the HII region (bottom 
panel) on the absorption feature. Here we adopt the same 
simplistic parameterisation as in Equation ([2]): an HII region 
of size Ri, surrounded by a homogeneously neutral IGM with 
neutral fraction xh- The dashed lines in the top panel rep- 
resent an IGM with xh = 0, and the solid lines represent 
the case with xh = 0.5 and Rb = 10 Mpc. IGM absorption 
produces a wider feature than that from a DLA. However, it 
becomes hard to visually discern absorption owing to a neu- 
tral IGM for N m > 10 21 cm -2 . This illustration suggests 
that it is difficult to measure xh from GRB050904, which 

3 To derive Equation J2]l, approximate a a (v) with a Lorentzian 
profile and ignore the redshift dependence that appears in the 
numerator of the integrand in Equation JTJ. We have checked 
that these approximations are accurate to better than 10% for 
relevant Rb- 



has N m » 10 216 cm" 2 ()Totani et al.ll2006T ). We return to 
this point in Section O 

DLAs with A r Hi > 10 19 cm -2 are associat ed with all 
but o ne GRB for which A^hi has been measured (|Chen et al.l 
l2007l V The cumulative distribution of Nm for the current 
sample of ps 30 GRB DLAs scales as Nm between 10 18 and 
10 21 ' 5 cm -2 , and about half of all GRBs for which A^hi ha s 
been measured have A^hi > 10 cm - 2 l|Chen et alj|2007h . 
However, it is not clear how Am should scale with redshift. 
It is plausible that, since the average galaxy becomes less 
massive with redshift, z > 6 galaxies should, on average, 
have weaker DLAs. For galaxies to reionise the Universe, 
the escape fraction of ionising photons must be appreciable, 
implying that sight-lines with Nm < 10 18 cm -2 must existQ 

The bottom panel in Figure [1] depicts the effect of bub- 
ble size on the GRB afterglow spectrum. The IGM absorp- 
tion for GRBs in bubbles with sizes between 1-10 Mpc is 
comparable. However, if the GRB sits in a large bubble with 
Rb ~ 50 Mpc, the absorption is strongly diminished. When 
xh = 0.2, roughly half of the skewers from GRBs in our 
simulations of reionisation sit in HII regions that are larger 
than 50 Mpc (Section ©. 



3 EFFECT OF PATCHY REIONISATION 

While the GRB itself does not ionise the IGM around it, 
earlier star formation from its host and neighbouring galax- 
ies can grow a large HII region. To model reionisation 

4 An optical depth of unity at the HI Lyman-limit requires A^hi = 
2 X 10 17 cm" 2 . 
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around a GRB, we employ three radiative transfer simu- 
lations that are each 186 Mpc on a side. These simulations 
are post-processed on a 1024 3 particle N-body field. Un- 
resolved halos are included with extended Press-Schechter 
merge r trees. These three sim ulations are described in de- 
tail in IMcQuinn et al] ((2 007b) [and the methodology is de- 
scribed in iMcQuinn et al. I i|2007al )], and they are meant to 
span the range of plausible morphologies for reionisation by 
galaxies. Here is a brief description of the models for the 
sources and sinks of ionising photons used in each of the 
three simulations: 

Model (i): All halos above the mass at which the gas can 
cool atomically (m coo i) contribute ionising photons at a rate 
that is proportional to their mass m. The scaling JVi on ~ m 
assumes that the massive SFR is proportional to the amount 
of gas within a galaxy. 

Model (ii): Halos more massive than m coo i contribute to 
the production of ionising photons, with the ionising lumi- 
nosity of the sources scaling as halo mass to the 5/3 power. 
This scaling is chosen to match the relationship between 
star formation efficiency and ga laxy mass that is obser ved 
in low-redshift dwarf galaxies (jKauffmann et all I2003T ) as 
well as the SFR fo und in theoretical studies that include su- 
pernova feedback (iDekel fc Wodl2003l ; lHernquist fc Springel 



120031 ; ISpringel fc Hernquistll2003l ) 

Model (iii): Absorption by minihalos shapes the morphol- 
ogy of reionisation and iVj on ~ m for m > m coo i . All mini- 
halos with m > 10 Mq absorb incident ionising photons 
out to their virial radii until they are photo-evaporated. We 
us e the fittin g form ula for the evaporation timescale given 
m llliev et all l|2005h . This timescale is roughly the sound- 
crossing time of a halo, or t cv — 100 (m/10 7 Mq) Myr 
l|Shapiro et al.ll2004 ). 

The normalisation of the function Ni on (m) in the three 
simulations is chosen such that reionisation is completed by 
z ~ 7. Given the uncertainties in / esc - the fraction of ionis- 
ing photons that escape and ionise the IGM - and in the SFR 
in high-redshift galaxies, there is a large range of possible 
normalisations. Fortunately, the morphology of reionisation 
when comparing at fixed xh dep ends only weakly on th e 
normalisation o f Nj on , a s sho wn in IMcQuinn et all (l2007al ). 
IMcQuinn et all (|2007ah and IMcQuinn et al.l (|2007rj| ) also 
showed that other effects such as thermal feedback on sub- 
Jeans mass galaxies, source duty cycle, and the number of 
recombinations have negligible impact on the morphology 
for reasonable models. In fact, studies have demonstrated 
that the morphology of reionisation is s haped principally 
by the c lustering of the ionising sources dFurlanetto et al.l 
l2004al fbl. 120051 ; IZahn et al.ll2007l ; IMcQuinn et al.ll2007ah . 

Figure [5] displays slices through simulations adopting 
reionisation models (i), (ii), and (iii). The white regions are 
ionised and the black are neutral. Model (ii) results in the 
largest HII regions because it has the most biased sources, 
whereas model (iii) produces the smallest bubbles, with the 
maximum bubble radius restricted to be roughly the mean 
free path for ionising photons to intersect a minihalo. 

The absorption profile of a GRB afterglow probes a line 
of sight through the IGM. Most of the absorption from the 
IGM occurs at < 70 Mpc from the host galaxy, and the neu- 
tral regions in the IGM that lie closer to the GRB contribute 
more absorption than those further away. It is clear from 
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Figure 3. Nine randomly selected 70 Mpc skewers originating 
at a GRB in the simulation volume. These sight-lines are chosen 
from a snapshot from the simulation of model (i) that has xh = 
0.5 and is from 2 = 7.3. 
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Figure 4. Transmission for the nine skewers displayed in Figure 
[3] Each transmission curve is calculated from the sight-line in 
Figure[3]that bears the same number. 



Figure [2] that different sight-lines encounter vastly different 
spatial distributions of xh- To further illustrate this point, 
Figure [3] displays nine randomly selected 70 Mpc skewers 
originating at a potential GRB site in the simulation volume. 
These skewers are from the z — 7.3 snapshot of the simula- 
tion using model (i) for which the volume-averaged neutral 
fraction is xh = 0.5 and the mass-averaged neutral fraction 
is 0.40 Henceforth, xh refers to the volume-averaged neu- 
tral fraction. To assign locations for GRBs in the simulation 



a The ionisation fronts are much narrower than our grid size for 
reionisation by POPII-like stars (as is assumed in these simu- 
lations). Resolution effects smear out the simulation ionisation 
field, making it deviate from a field of 0s and Is. This smearing 
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volume, we assume that the -/Vi on of a halo is proportional 
to the GRB rate (which is true if the GRB rate traces the 
massive SFR and / csc is independent of halo mass). 

Figure [3] plots the IGM transmission redward of Lya 
for the nine GRB sight-lines shown in Figure [3] To compute 
these curves, we use Equation (1) and the density, ionisation, 
and velocity fields from the simulation. The observed trans- 
mission across these six lines of sight varies enormously, and 
this variance needs to be accounted for in measurements of 
xh using GRBs. These calculations do not account for the 
evolution in xh along the line of sight owing to light-travel 
effects, but instead compute the absorption from a skewer 
using a snapshot that is fixed in time. In the time for light to 
travel 50 Mpc, the value of xh in the simulations of model 
(i) changes by 0.06, 0.08, and 0.13 when xh = 0.7, 0.5, and 
0.2, respectively. A proper treatment of the evolution of xh 
along a sight-line would hardly affect our conclusions. 

For the nine sight-lines in Figure^ the best fit (xh, Rb) 
to the afterglow transmission redward of Lya are respec- 
tively (0.31,6.6), (0.35,14), (0.47,2.3), (0.51,0), (0.28,0), 
(0.21, 10), (0.32, 22), (0.45, 10), and (0.43, 2.5)0 While these 
two parameters do not fully characterise the distribution of 
xh, the actual value for R b of the GRB host bubble corre- 
sponds roughly to the size given by these fits and the fitted 
value for xh is within 0.3 of the global value. In reality, 
a DLA will also affect the damping wing profile, and its 
absorption must be accounted for simultaneously. We incor- 
porate this extra complication in the following section. 

The top left panel in Figure [5] plots the probability dis- 
tribution of GRB host HII bubble radii, P(R b ). The GRB 
host bubble radius is defined as the distance it requires for 
a skewer to cross three simulation grid cells that are more 
than 1% neutral. The dashed, solid, and dot-dashed curves 
are Rb P(Rb) for xh = 0.2, 0.5 and 0.7, respectively. The 
thick curves are for model (i), and the medium width and 
thin solid curves (shown for only xh ~ 0.5) adopt model (ii) 
and model (iii), respectively. These curves demonstrate that 
the function P(Rb) depends more strongly on xh than on 
the reionisation model. When xh ~ 0.5, the typical GRB 
resides in a bubble of size Rb ~ 10 Mpc, and when xh ~ 0.2 
this number increases to Rb ~ 60 Mpc. 

The top right panel in Figure [5] shows the distribution 
of xh along the line of sight to a GRB at fixed global xh- 
To create this PDF, we generate 10 6 skewers that originate 
at GRBs in the simulation volume and then tabulate xh 
along each skewer, weighting points as 1/r 2 where r is the 
distance from the source - roughly the same weighting that 
is relevant to the damping wing optical depth at observed 
frequency u a /(l + z g ) (eqn. [2}. We do not include points 
inside the host bubble in this tabulation. 

The wide breadth of the xh distributions indicate that 
a single GRB cannot be used to directly measure the global 
value of xh to a precision of Sxh ~ 0.3. Despite the scatter 
in xh, a single GRB can be used to detect whether the IGM 



does not appreciably affect the results of the calculations reported 
here. 

6 Note that the two skewers which are best fit with Rb = are 
anomalous, and the total fraction of lines of sight that are best 
fit with Rb = is a much smaller fraction than suggested by this 
sample of nine sight-lines. 
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Figure 5. Dashed, solid, and dot-dashed curves represent xh = 
0.2, 0.5 and 0.7, respectively. The thick curves are calculated us- 
ing model (i), and, in the top panels, the medium width and thin 
solid curves (shown only for xh = 0.5) use model (ii) and model 
(iii). Top Left Panel: The PDF of R b constructed from 10 6 sight- 
lines, in which Rb is defined as the length of a skewer originating 
at a GRB that passes through three grid cells with xh > 0.01. 
Top Right Panel: The PDF of xh measured directly from 10 6 
lines of sight as discussed in the text. Bottom Panels: The same 
as the top panels except that the PDFs of xh and Rb are con- 
structed from fitting the afterglow transmission curves computed 
from 1000 different simulation sight-lines with these parameters. 



is neutral. We focus on the probability of a GRB allowing 
the detection of reionisation in the next section. 

It is possible that GRBs occur in only the most massive 
galaxies (rather than in a typical one, as we have assumed 
thus far). If only halos with m > 10 10 Mq are able to pro- 
duce GRBs with a rate proportional to the halo mass, again 
weighting these halos by their SFR, we find that the average 
bubble size increases by a comparable amount to the size in- 
crease between model (i) and model (ii) shown in Figure [5] 
However, we find that the low-end tail of P(Rb) starts to 
disappear with increasing m, even though the peak of the 
PDF does not shift significantly. 

Thus far, we have discussed the distribution of xh and 
Rb measured directly from the simulations. For a GRB after- 
glow, xh and Rb are measured by fitting these parameters 
to the afterglow damping wing (even though this parameter- 
isation for the distribution of xh along a sight-line is highly 
simplified) . The bottom panels in Figure [5] show the PDF 
of Rb and xh from fits to 1000 transmission curves com- 
puted from GRBs in the simulation of model (i) . We ignore 
the effect of DLA absorption in these fits. These PDFs are 
notably different from those in the top panels - favouring 
smaller Rb and lower xh- While this procedure returns bi- 
ased values for Rb and xh relative to the values we estimate 
in the top panel, we find that this two parameter model pro- 
vides an excellent fit to the absorption profile (as we discuss 
in Section [4]). 



6 M. McQuinn et al. 




8900 



9000 9100 

X (Angstroms) 



Figure 6. Afterglow spectrum of a hypothetical 



9200 



6.3 burst, 



computed with xh = 
curves), Nm = 10 20 



0.5 and TVjji = 10 19 cm" 2 (thick solid 
cm -2 (medium-width solid curves), and 
Nm = 10 21 cm -2 (thin solid curves) for the specified Ri,. The 
dashed curves are the best fits for a model that includes only 
DLA absorption, fitted to AA = 200 (1 + z g ) = 1500A. The nor- 
malisation of the two sets of curves is adjusted to aid viewing. 



4 FITTING THE DAMPING WING IN GRBS 

This section quantifies the circumstances under which a 
GRB afterglow can be used to detect a neutral IGM. Previ- 
ous sections have illustrated how the patchiness of the IGM 
and the strength of the host DLA complicate this measure- 
ment. In addition, there are other uncertainties such as the 
power-law slope and amplitude of the intrinsic GRB spec- 
trum, the precise redshift of the GRB and of the DLA, and 
dust absorption and metal line contamination from the host 
galaxy or from intervening syste ms. A systematic d iscus- 
sion of these issues is presented in iTotani et al.l (|2006l ) . To 
rule out an ionised IGM with a GRB afterglow, a multi- 
parameter fit to its spectrum that assumes an ionised IGM 
and that accounts for the aforementioned uncertainties must 
provide a poor fit. 



4.1 Fits to a Toy Absorption Model 

A toy example of such a multi-parameter fit that accounts 
for some of these uncertainties is shown in Figure [6] In this 
figure, the simple case of a bubble with size Ri, in a ho- 
mogeneous IGM with neutral fraction xh is shown. The 
solid curves are absorption models that have the intrin- 
sic parameters f3 = 1.25, z g = zdla = 6.3Q xh = 0.5 
and AT H i = 10 19 cm" 2 (thick curves), Nm = 10 20 cm"' 
(medium-width curves), or A^hi = 10 21 cm -2 (th 



)H 



nn curves 



7 The assumption that z g = zdla is probably reasonable for 
GRBs in high-redshift galaxies, wh ere the circular velocit y of 
these halos is typically ~ 25 km s — 1 l lBarkana & Locb 2004TI, 

8 We include the prior 0.75 < (3 < 1.75 in the fits: /3 will likely be 
constrained by earlier photometric observat ions redward of th e 
absorption, as was the case for GRB050904 llTotani et alj|2006h . 
In addition, we exclude wavelengths within 400 kms" 1 of the 
host galaxy because such regions could be affected by Lyo emis- 



We then fit to these model curves. Our fitting weights all fre- 
quency bins equally. The presence of night-sky OH lines and 
the wavelength dependence of the CCD sensitivity should 
make the weighting slightly non-uniform in a real observa- 
tion. 

The dashed curves in Figure [S] are fits over AA = 
200 (1 + z g ) to a model that only include DLA absorption. 
This model fits for the parameters Nm, /3, and A. The ab- 
sorption owing to a substantially neutral IGM is more easily 
fit with just a DLA model as Rt, increases or as A^hi in- 
creases. If we had fit over a larger AA, the dashed and solid 
curves would differ by less at wavelengths where the absorp- 
tion is important. The fits in Figure [6] assume that zdla is 
known from metal lines, as was the case for GRB050904. For 
GRBs in more metal-poor galaxies or for GRBs with smaller 
DLAs, this may not always be possible. When fits include 
zdla as a free parameter, a neu tral IGM is much harde r to 
distinguish from an ionised one (jBarkana fc Loebll2004D . 

To quantify the capability of distinguishing different ab- 
sorption models, we define the measure 



(Ax 2 ) = E 



[F2(AQ-Fi(AQ] 2 
<7(A,) 2 



(3) 



where Fi and F2 are the afterglow fluxes for models 1 and 
2. This quantity represents the average difference in \ 2 be- 
tween two models for a spectrum that has noise cr(Aj) in 
channel Ai and that has N wavelength channels, assuming 
the data is drawn from model 1 or 2. 

When comparing a model with a neutral IGM to an 
ionised one, two models can be distinguished at X confi- 
dence level (C.L.) if X of the total likelihood is contained 
between the maximum likelihood in the neutral IGM model 
and exp[— (Ax 2 ) /2] of that likelihood value (assuming Gaus- 
sianity). For a 4 parameter fit, two models that differ by 
(Ax 2 ) = 8 signifies that model 1 is, in the mean, pre- 
ferred at 91% C.L. over model 2 ((Ax 2 ) = 5, 70% C.L. 
and (Ax 2 ) = 12, 98% C.L.). The significance levels at fixed 
Ax 2 for fits with three and five parameters are similar. In 
what follows, we quote (Ax 2 ) in terms of AA; = 3 A and 
err — 0.1, where AAi is the width of a frequency channel 
and ajr is the standard deviation in the flux in each spectral 
channel in units of A. These are approximately the values 
of AAi and ajr for the Subaru FOCAS spectroscopic obser- 
vation of GRB050904. 

First, we consider the case in which the absorption spec- 
trum is parameterised by Nm ~ 10 19 cm -2 , xh = 0.5, 
and R b = 10 Mpc (similar to the fits in Fig. gj). 8 We fit 
to AA = 50 (1 + z g ) A redward of source-frame Lya. This fit 
results in (Ax 2 ) = 5 x (3A/AAj) x (0.1/<7;r) 2 between the 
model for a neutral IGM and the model for an ionised IGM. 
Next, we fit to a model with the same specifications as the 
previous fit but with Nm ~ 10 20 cm" 2 [Nm ~ 10 21 cm" 2 ]. 
This fit results in (A X 2 ) = 3 [0.2] x (3A/AA*) x (0.1/W) 2 . 
Therefore, AAi < 3 i and ajr < 0.1 are required for any 
hope of detecting a neutral IGM for the Nm ~ 10 19 cm" 2 , 
xh ~ 0.5, and = 10 Mpc case, and even better sensitivity 
is required for the other two examples. 

Let us investigate how the quoted (Ax 2 ) depend on the 



sion from the galaxy as well as resonant absorption by infalling 
material. 



Probing the Neutral Fraction of the IGM with GRBs during the Epoch of Reionisation 7 



assumptions we have made. For this discussion, we again 
assume an absorption model with xh = 0.5, Rb ~ 10 Mpc, 
and iVm ~ 10 20 cm -2 , and we assume a measurement with 
A\ fa 3 A and a? ~ 0.1. 

• If we restrict /3 to 1.0 < /3 < 1.5 (rather than to 0.75 < 
P < 1.75, as before), (A X 2 ) increases from 3 to 4. If we fix 
P = 1.25, then (Ax 2 ) = 5. 

• If we allow zdla to vary and fix (3 = 1.25, (Ax 2 ) de- 
creases to 1. 

• If we instead fit the first A A = 200 (1 + z g ) A redward 



of Lya rather than AA 
from 3 to 6. 



50 (1 + z g ) A, (Ax ) increases 



The increase in (Ax 2 ) with bandwidth owes to tighter con- 
straints on A and (3, the parameters that determine the spec- 
trum far redward from the Lya absorption feature. These 
tighter constraints breaks degeneracies between these pa- 
rameters and the parameters that determine the absorp- 
tion. The appropriate choice of bandwidth will not always 
be clear-cut; fits should not include wavelengths where the 
dust absorption, the instrumental response, and other fac- 
tors are not understood to better than the difference be- 
tween the two absorption models. Already for AA « 300 A, 
such calibration is necessary to the 1% level (Fig. |6}. In 
addition, even uncertainties in the cosmological parameters 
make a difference in distin guishing models at the percent 
level (iMiralda-Escudei 1 19981 ) . 
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Figure 7. PDFs of (Ax 2 ), where (Ax 2 ) is the average difference 
in x 2 between the "true" model (which fits xh and Rb) and a 
model with only DLA absorption. These PDFs are generated by 
fitting 1000 sight-lines from the simulation of model (i) at the 
specified xh- Note that the parameter y that is plotted on the 
abscissa is scaled to be independent of S/N and R. The left panel 
shows the dependence of the PDF on JVhi, using spectra from a 
snapshot with xh = 0.5, and the thick curves in the right panel 
show the dependence on xh, assuming TVjii = 10 20 cm -2 . The 
thin black curve in the right panel folds in the observed distribu- 
tion of ./Vhi and uses spectra with xh = 0.5. A neutral IGM is 
detected on average at > 84% confidence level when (Ax 2 ) > 8 
(marked by the vertical dashed lines for R = 3000 and S/N = 10). 



4.2 Fits to the Simulated Absorption Spectra 

Thus far, we have considered fits to a toy model for IGM 
absorption. To realistically model IGM absorption we cal- 
culate mock spectra from simulation skewers of reionisation 
model (i) and we set f3 = 1.25. We fit this mock data with 
different absorption models in what follows. 

Despite the complicated distribution of neutral hydro- 
gen along sightlines, a simple model with just the parameters 
Rb and xh (in addition to parameters that account for DLA 
absorption and the intrinsic afterglow spectrum) typically 
provides an excellent fit to the mock data. This is surprising 
because of the highly simplified nature of this model. The 
(Ax 2 ) between the true model - the absorption profile cal- 
culated from the simulation - and the best fit neutral IGM 
model, which fits for Rb and xh, typically range between 
10" 3 - 10" 7 x (3A/AA0 x (0.1/cr^) 2 for N m = 10 20 cm" 2 . 
These small values mean that ajr ~ 10~ 2 — 10 -4 is required 
to rule out this 2-parameter IGM absorption model (mea- 
surements ~ 10 — 1000 times more sensitive than those for 
GRB050904). 

Figure [7] shows the distribution of (Ax 2 ) between fits 
to the simulated spectra with the parameters Rb, xh, Ahi, 
A and /3 versus fits with just the last three parameters (i.e., 
a DLA and an ionised IGM). These fits are to a wave- 
length interval that spans 50 (1 + z) A redward of rest- 
frame Lya. If the fits with xh = are disfavoured by 
the data (i.e., if (Ax 2 ) > 8), then a neutral IGM will 
be favoured by the data. The parameter y that is plotted 
on the abscissa in Figure [7] is scaled to be independent of 
S/N (we define S/N = Up 1 ) and spectrograph resolution R 
(R = Ai/AAi). For an observation like the FOCAS obser- 
vation of GRB050904 (R = 3000 and S/N = 10), a DLA 
with iVm < 10 20 cm -2 is needed to have a good chance of 



detecting a neutral IGM when xh ~ 0.5 (left panel in Fig. 
UJ. For xh < 0.5, a higher signal-to-noise measurement is 
required. 

The results in Figure [7] depend on the parameters we fit 
as well as the wavelength range included in the fit. If we fix 
the parameter j3 = 1.25 instead of fitting for it, the PDFs 
in this figure are, for the most part, unchanged. If we fit a 
larger range in wavelength of A A = 200 (1 + z) A, the fits 
improve slightly - shifting the histogram by factor of ~ 1.5. 

The thin solid curve in the right panel of Figure [7] is the 
PDF of (Ax 2 ) for spectra computed for xh = 0.5 and where 
we have folded in the distribution of DLA Nm values ob- 
served in lower redshift GRB-host systems. We assume that 
the cumulative distribution of DLAs is a power-law of the 
form iVSi 3 for N m < 10 21 5 cm , which is consistent with 
observations (|Chen et all 120071 ). Half of GRB DLAs have 
Am > 10 21 ' 5 cm~ 2 , but we do not include these bursts in 
this computation. These bursts will contribute to the small 
(Ax 2 ) tail of the PDF. 

Let us assume there exists an observation with similar 
sensitivity to the Subaru FOCAS spectrum of GRB050904, 
that the distribution of DLAs is the same as found at lower 
redshift, that the redshift of the GRB is known via metal 
lines, and that the GRB originates from a redshift at which 
xh ~ 0.5. The spectrum of this GRB can be used to de- 
tect a partly neutral IGM at 97% C.L. ((Ax 2 ) > 12) ap- 
proximately 5% of the time (and, for an observation with 
3 times the sensitivity, « 25% of the time) if the fit uses 
AA = 50 (1 + z) 10 If the fit uses AA = 200 (1 + z) A, this 

9 These numbers assume that bursts with Nm > 10 21 ' 5 cm -2 
do not provide a positive detection of a neutral IGM. 
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becomes w 10% of the time (and, for an observation with 3 
times the sensitivity, m 30% of the time) . 

An earlier spectroscopic pointing could result in higher 
S/N values than the fiducial value of 10 because the optical 
afterglow fades as a p ower-law in time with slope » —1.2 
i|Liang fc Zhanell2006h . The same Subaru observation taken 
» 8 hr after GRB050904 would have resulted in 10 x the 
sensitivity and a detection of reionisation the majority of 
the time. In 2013, JWST is projected to be in orbit. Let us 
assume a z = 8 GRB is observed with a 10 4 s integration us- 
ing the JWST spectrograph NIRSpec operating in its high- 
est resolution mode of R — 2700. To detect the GRB with 
S/N = 10 in each frequency channel requires a minimum 
flux of F(X) = 2 x 10~ 19 erg cm -2 s _1 A _1 This mini- 
mum flux is approximately five times the sensitivity of the 
Subaru FOCAS observation of GRB050904. However, while 
ground based observations can potentially target GRB af- 
terglows minutes after their trigger, JWST will only be able 
to slew to high-redsh ift afterglows ~ 1 day after the burst 
l|Gardner et alll2006j ). 

Even though the calculations in this section were for 
z » 6 — 8, similar conclusions hold for higher redshifts reion- 
isation scenarios. The amount of absorption does increase 
with increasing z g because the damping wing optical depth 
scales as tn (1 + z g ) 3 ^ 2 (eqn. [2}. However, the morphology 
of reionisation depends weakly on redshift for stellar reioni- 
sation scenarios - bubbles have a similar comoving size dis- 
tribution at a fixed xh indepen dent of when reionisation 
happens (McQ uinn et al. 2007a). For reasonable redshifts 
over which reionisation could occur, the amount of damping 
wing absorption predicted by our simulations is similar and 
our conclusions do not change. 

As a final point, the distribution of the best-fit Rb and 
xh from the fits that include DLA absorption are very sim- 
ilar to the distribution in Figure 5 in which the effect of a 
DLA was ignored. Therefore, the presence of a DLA does 
not significantly bias the best fit values for Ri, and xh- 

5 GRB050904 

The z = 6.3 GRB, G RB050904, is the burst with the high- 
est identified redsh ift l|Kawai et al.ll2006l ; lHaislip et al.ll2006l ; 
iTotani et al.l 120061 ) . A 4 hour observation of this burst was 
taken using the FOCAS spectrograph on Subaru 3.4 days af- 
ter the prompt gamma ray emiss ion. The spec t rum f rom this 
observation is shown in Figure 151 ITotani et"afl (2006) showed 
that this spectrum is well-fitted with a DLA of column den- 
sity N m « 10 21 ' 6 cm -2 and that this GRB disfavours addi- 
tional absorption owing to a neutral IGM, deriving the con- 
straint xh < 0.6 at 95% C.L. This section investigates the 
constraints on xh from GRB050904 if a patchy reionisation 
process consistent with simulations is taken in to account. 

First, we have re-performed the analysis of lTotani et al.l 

(2006), accounting for galactic absorption of E(B — V) = 
0.06 and ma s king t he same regions of the spectrum as in 
ITotani et ail (2006) because of metal line contamination. 
(The thick horizontal lines at the top of Figure[8jindicate the 
wavelength ranges that were fit.) The solid absorption curve 

10 http:/ /www. stsci.edu/jwst/instruments/nirspec/ 
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Figure 8. Binned afterglow spectrum of GRB050904 (each 
2.67 A wavelength channel has been binned into 8 A pixels), as 
well as the best fit models for xh = (solid curve) and xh = 1 
(dashed curve). The bottom curves are the 1-cr errors on each 
binned pixel with an offset of —5 X 10 — 19 , and the thick hori- 
zontal lines at the top of the figure indicate the wavelengths that 
were included in the fits. Both fits require a DLA with column 
density N m s» 10 215 cm" 2 . 

in Figure [8] is the best fit for an ionised IGM (xh = 0) with 
a DLA and the dashed curve is the best fit with a DLA plus 
a fully neutral IGM. Both fi ts to the spec t rum e mploy the 
same estimated errors as in ITotani et ail (2006) and min- 
imise x 2 with two parameters, A^hi and A, holding all other 
pa rameters fixe d . In p articular, we fix (5 = 1.25 as was done 
in ITotani et alj (|2006l ). This value for (3 is consistent with 
(3 = 1.25L0.25, which was derived in iTagliaferri et all (|2005l) 
and is consistent with the value derived in lKann et alj (|2007J) 
from all publicly available data(3 The value of A\ 2 is not 
significantly changed if we instead fit with j3 adopting the 
prior 1.0 < /3 < 1.5 (the best fit for both models then prefers 
(3 = Vj, and Ax 2 shrinks significantly if no prior is placed 
on/30 ' 

The \ 2 values for the two models displayed in Figure 
[8] are 277.3 for the ionised model and 287.0 for the neutral 
case0 Even though these two 3-parameter models are dif- 

11 The ITagliaferri et alj <2005t) value for (3 was derived from pho- 
tometric observations with the Very Large Telesco pe taken 1.2 
days after the prompt emission. Totani et al. (2006) argued that 
the softness of the derived in ITagliaferri et al.l ||2005|) indicates 
that (3 was determined by the power-law index of the electrons 
and, therefore, should no t evolve significantly b etween the time 
of the measurement from ITagliaferri et alj | |2005|) and their mea- 
surement (i.e., other breaks in the synchrotron spectrum are at 
longer wavelengths). 

12 If we ignore photometric constraints on (3 and just fit it as a 
free parameter, the ionised IGM model prefers (3 = —0.6 at 98% 
C.L. (Ax 2 = 10) over a model with = 1.25. While it is hard 
to imagine (3 evolving to this small value from j3 = 1.25, this 
nevertheless may indicate that there is a problem with the model 
for the data. Fits to the spectrum that do not place a prior on (3 
do not prefer ionised IGM m odels over neutral o nes. 

13 To emulate the results of lTotani et al. I J2006h . the error bars 
have been renormalised by the factor of 0.848 from their origi- 
nally estimated value such that the x 2 /d.o.f. 1 for the best fit 
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ficult to distinguish visually, it is the case that Ax 2 = 9.7, 
which at face value implies that the ionised model is pre- 
ferred over the neutral model at 98% C.L. (jTotani et al,l 
2006). Interestingly, much of the contribution to A\ 2 be- 
tween the neutral and ionised models comes from wave- 
lengths that are not significantly affected by IGM and DLA 
absorption. If we tabulate Ax 2 for the best fit models in 
Figure [8] only including the wavelengths within 200 A of 
Lya - the wavelengths most affected by absorption - in the 
summation, the two model fits differ by Ax 2 = 4. If we in- 
clude 700 A redward of Lya, the two model fits in Figure [8] 
differ by A^ 2 = 5. The significant contribution to A^ 2 from 
redward of the Lya absorption stems from the fact that the 
neutral IGM model prefers a 2.4% higher normalization to 
compensate for the additional effect of IGM absorptions. As 
a consequence, the unabsorbed part of the spectrum needs 
to be modelled to a relative precision of a percent over a 
wavelength range with A A/A m 0.1 for a correct interpreta- 
tion of the data. 

To distinguish models that differ at the 2% level over 
hundreds of angstroms posits that intergalactic metal and 
sky lines have been properly masked, that the dust redden- 
ing of the Milky Way and the host galaxy has been accu- 
rately accounted for, that the relative calibration of and er- 
rors on the measured flux are accur ate, and that t he m odel 
for the intrinsic spectrum is correct. iTotani et al.l l|2006l ) in- 
vestigated many of these uncertainties and argued that the 
FOCAS spectrum of GRB050904 can be used to distinguish 
the models in Figure ISP^I 

The above analysis did not fit for Rb and xh, which 
must be done to place constraints on xh for GRBs dur- 
ing the Epoch of Reionisation. However, while GRB050904 
may favour an ionised universe over a neutral universe, it 
is more difficult to place a constraint on intermediate xh 
with GRB050904, especially in light of the large bubble sizes 
during reionisation. If we fit a model that fixes xh = 0-5 
(abou t as low an xh as is plausible at z = 6.3 ijLidz et al.l 
120071 ')') and R b = 10 Mpc [30 Mpc] - again fitting with A 
and A^hi and fixing f3 — , then the Ax 2 between a model 
with an ionised universe and this model is 3.1 [2.5], which is 
not statistically significant. To place any constraint on xh 
for GRB050904, one must account for the covariance of Rb 
and xh (i.e., lower xh have smaller Rb). We have not taken 
the additional step of accounting for this covariance in our 



analysis because it will not lead to interesting constraints on 
x H from GRB050904. For future bursts that prefer no IGM 
absorption, such an anlysis will be essential. In any case, ac- 
counting properly for patchy reionisation would significantly 
weaken the 95% C .L. constraint xh < 0.6 at z — 6.3 that 
ITotani et all ({2006) derives from GRB050904. 



6 THE DAMPING WING IN QSOS AND 
GALAXIES 

GRBs are not the only beacons in which the signature of 
a neutral IGM can be observed in their continuum emis- 
sion. Rather than wait for a high-redshift GRB to occur, 
hundreds of z > 6 galaxies and QSOs have already been 
found. Perhaps these objects can be used to detect a neu- 
tral IGM. However, since the observed population of QSOs 
and galaxies are more biased tracers of the high-redshift 
Universe than are GRBs (assuming that GRBs trace star 
formation), they sit in larger bubbles such that the damp- 
ing wing absorption is smaller, on average. The effective size 
of an HII region for one of the known z ~ 6 quasars, assum- 
ing that the IGM is significantly neutral at this redshift, is 
predicted to be > 50 Mpc. This number accounts for these 
rare objects being in the most overdense regions in which 
reionisation occurs earlier (|Lidz et al.ll2007l ). Such large HII 
regions make searches that target the red damping wing in 
the highest redshift quasars hopeless PI 

Galaxies will sit in smaller HII regions than QSOs. The 
average bubble size around galaxies is larger than those for 
GRBs. As noted in Sj3] the difference in the typical bubble 
size is typically less than a factor of two between halos of 
m ~ 10 10 " 11 Mq (roughly the mass of halos of spectroscop- 
ically confirmed z > 6 galaxies) and GRB-hosting halos. 

However, there are two significant complications. First, 
because galaxies are much fainter, the continuum emission 
for high-redshift galaxies is difficult to detect spectroscopi- 
cally with existing facilities. Stacking galaxy spectra is nec- 
essary to have any hope of detecting a damping wing feature 
with current data. Another significant complication is that 
the continuum of galaxies is far from a simple power-law. 
Such uncertainties must be accounted for in any analysis 
that attempts to obtain xh from the damping wing absorp- 
tion profile of a galaxy. 



ionised model, where the number of degrees of freedom (d.o.f.) 
equals 268 — 3. ITotani et al. I <2006h uses this rescaling to correct 
for uncertainty in the errors. The residuals for the best fit ionised 
model are close to a Gaussian with th e correct width after this 
pr ocedure (Fig 3 inlTotani et al.ll2006T) . 

Tot ani et al. also used the transmission in the Ly/3 for- 

est as an additional handle to discriminate between the best fit 
for a model that assumed DLA absorption and an ionised IGM 
and the best fit model that assumed the absorption owed entirely 
to a neutral IGM. The later fit required the redshift of the host 
galaxy to be significantly redward of the best fit redshift of the 
DLA absorption lines, which was inconsistent with the transmis- 
sion in the Ly/3 forest. This example is a rather exceptional case, 
and for most bursts Ly/3 will not provide an additional handle. 
In addition, patchy reionisation all ows for t ransmission blueward 
of GRB-frame Ly/3. Therefore, the lTotani et al" I <2006ft argument 
can be applied only when the model value for xh is larger than 
the fraction of pixels that have transmission. 



7 CONCLUSIONS 

GRBs are the most luminous sources at high redshifts, and 
their smooth power-law spectra are ideal for isolating the 
effects of absorption owing to a neutral IGM. However, ob- 
servations must separate the impact of IGM absorption from 
that of a DLA to detect a neutral IGM. If no damping wing 
feature from IGM absorption is detected in the spectrum of 



15 Rather than analyse wavelengths redward of the QSO 
Lya line, damping wing absorption from neutral patches in 
the IGM would also affe ct the Lya forest of these QSOs 
iMesin gcr & Haiman 2007). Owing to little forest tr ansmission 
at z > 6, this feature would be difficult to detect llLidz et al.l 
l2007t) . 
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a high-redshift GRB, one must be careful to conclude that 
reionisation is complete at the redshift of interest. 

We have shown that there is a wide probability distri- 
bution of HII region sizes, and that there is large variation 
in the distribution of xh along different sight-lines. A non- 
detection of neutral hydrogen from a GRB afterglow spec- 
trum might arise because the GRB host galaxy sits within 
a large HII region. If absorption owing to a neutral IGM 
is detected, it will be impossible to infer xh from a single 
GRB to better than Sx h ~ 0.3 because of the patchiness of 
reionisation. 

Assuming an observation with similar sensitivity to the 
Subaru FOCAS spectrum of GRB050904, that the distribu- 
tion of DLAs is the same as found at lower redshifts, and 
that the redshift of the GRB is known, a GRB from a red- 
shift at which xh ~ 0.5 can be used to detect a partly 
neutral IGM at 98% C.L. « 10% of the time (and, for an 
observation with 3 times the sensitivity, ps 30% of the time). 
If xh < 0.5, these percentiles for detection are even smaller. 
Weaker DLAs enhance the probability of detecting a neutral 
IGM, but too weak of a DLA may prevent a precise redshift 
determination, which is essential for tight constraints on xh- 

Since the z = 6.3 burst GRB050904 has a DLA with 
Ahi = 10 21 ' 6 cm -2 , the absorp tion on the red sid e of the 
line is dominated by the DLA (|Totani et al.ll2006l ). While 
this burst may fav our a model with a n ionised universe over 
a neutral universe l|Totani et al.ll2006h , a weaker DLA is nec- 
essary to be able to constrain xh in the spectrum of a high- 
redshift GRB. This is particularly true if the GRB occurs 
within a large HII region. 

GRB050904 was observed spectroscopically by the Sub- 
aru telescope 3.4 days after the prompt gamma ray emission. 
If this afterglow had been observed hours after the burst, 
the flux would have been more than an order of magni- 
tude larger. To detect a neutral IGM, it is crucial for optical 
and near-infrared spectrographs to observe candidate high- 
redshift GRBs as soon as possible after the prompt gamma 
ray emission. A high signal-to-noise ratio is critical to distin- 
guish IGM absorption from that arising from a DLA. Such 
an observing programme is worthwhile given the promise 
that GRBs have as probes of the epoch of reionisation. 
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